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CREEP  BEHAVIOR  OF  UTAH  OIL  SHALE 
SUBJECT  TO  UNIAXIAL  LOADING3 

By 

Ken  P.  Chong1 ,  George  F.  Dana2  and  Jui-Lin  Chen 


ABSTRACT 

This  paper  presents  results  of  a  study  on  the  creep  behavior  of 
Utah  oil  shale.   A  Conbel  Model  355  pneumatic  driven  testing  machine 
is  used.   The  set  of  duplicate  test  specimens  required  for  creep  test- 
ing were  cut  from  a  single  oil  shale  layer  using  a  wire  saw  to  avoid 
any  surface  damage.   A  rheological  model  was  developed  for  creep  behav- 
ior of  oil  shale  as  a  function  of  stress  level  and  organic  content. 
Data  from  creep  testing  and  Fischer  assay  analyses  were  used  to  demon- 
strate correlation  between  various  stress  levels  and  organic  contents 
for  samples  taken  from  the  Mahogany  Zone  of  the  Parachute  Creek  Member 
in  Utah's  Cowboy  Canyon. 


INTRODUCTION 
Oil  shale  resources  in  the  United  States  are  important  as  an  al- 
ternative energy  source.   Although  a  number  of  development  procedures 
are  being  investigated,  the  need  to  modify  the  rock  during  mining  or 
in-place  is  common  to  all  of  them.   Knowledge  of  the  mechanical  prop- 
erties of  oil  shale  is  essential  to  any  form  of  development,  including 
room-and-pillar  mine  design.   One  of  these  mechanical  properties  of  oil 
shale  is  its  creep  behavior. 

In  this  paper,  uniaxial  creep  tests  were  performed  to  investigate 
the  effects  of  different  stress  levels  and  organic  contents  on  the 
creep  behavior  of  oil  shale.   Each  set  of  test  specimens  was  prepared 
from  a  single  layer  of  oil  shale.   The  lateral  homogeneity  of  oil  shale 
(7)  insures  that  the  sets  were  duplicate  samples.   All  specimens  were 
cut  perpendicular  to  the  bedding  planes.   Statistical  analysis  indicates 
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that  the  creep  properties  of  oil  shale  are  strongly  dependent  on  organic 
contents  and  various  stress  levels. 


THEORETICAL  CONSIDERATIONS 

The  objective  of  this  section  is  to  provide  the  theoretical  back- 
ground and  the  type  of  rheological  model  (1,  5,  6,  8,  9)  used.   Oil 
shale  is  nearly  linear  elastic  under  small  strains  as  measurements  of 
load-displacement  relationship  reveal  (2,  3,  4).   However  under  constant 
loading  it  creeps  slightly.   Thus  the  model  that  was  chosen  is  for  the 
materials  which  retains  initial  linearity  between  load  and  deflection 
and  a  non-linear  relationship  depending  on  time. 

Based  on  the  observation  of  the  raw  data  on  e  vs  time  plots  (Fig.  1) 
the  following  nonlinear  model  (8)  is  used: 

e  (0  ,  S,  t)  =  e°(0  ,  S)  +  e+  (0  ,  S)  tn (1) 

c  c  c 

where,  e  =  total  strain  perpendicular  to  bedding  planes  =  £33* 

e°  =  initial  strain; 

e  =  creep  related  strain; 

t  =  non-dimensionalized  time; 

0  =  organic  content  expressed  as  a  volume  fraction; 

S  =  stress  level  expressed  as  a  fraction  of  ultimate  stress;  and 

n  =  coefficient  depending  on  0  and  S. 

To  determine  the  functions  e°,  e   and  n,  the  following  procedure 
is  used: 

Taking  the  logarithms  of  Eq.  (1)  yields 

log  (e  -  e°)  =  log  e  +  n  log  t (2) 

Fig.  2  shows  a  typical  plot  of  log  (e  -  e°)  versus  log   t   result- 
ing in  a  straight  line  with  slope  n  and  the  intercept  at  unit  time  is 
c+ .      On  the  other  hand,  e°  is  the  intercept  at  0   time  taken  from 
Fig.  1. 

Based  on  the  present  uniaxial  testing,  a  nonlinear  pseudo  3-D  creep 
model  can  be  constructed,  relating  the  strains  to  stresses. 

For  pure  axial  compression,  the  stress  tensor  is 
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Fig.    1.      Longitudinal   Strain  vs.    Time. 
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For  simplicity,  let  o ?   =  o 


(4) 


where  axis  1,  2  indicate  directions  along  the  bedding  planes  and  axis  3 
is  perpendicular  to  the  bedding  planes,  forming  a  right-handed  Cartesian 
coordinate  system.   The  nonlinear  strains  are  (8) 

e33  (0  ,  S,  t)  =  Fi  a  +   F2  a2  +  F3  a3 

£22  (0c,  S,  t)  =  en(Oc,  S,  t)  =  v31  e33   (0c»  S,  t) 

=  (Fi  -  Gi)a  +  (F2  -  G3)a2+  (F3  -  G^o* 


^32  -  £31  -  ^21-0  

From  previous  investigations  (4) , 


(5) 


v3i  =  -0.2747  +  1.162  0  +  0.642  S 
dl  c 


(6) 


where  Fi ,  F2  ,  F3  ,  Gi  ,  G3,  Gi+  are  kernel  functions  involving  the  var- 
iables 0  and   t.   These  6  functions  can  be  solved  from  the  following 
6  simultaneous  equations: 


e33  (0  ,  0.25,  t)  =  Fi  o  +  F2  a2  +  F3  c3  .  . 

C                           a.  cd  3. 

e33  (0c,  0.50,  t)  =  F1   ab  +  F2  a2  +  F3  a3  .  . 

e33  (0c,  0.75,  t)  =  Fi  ac  +  F2  a2  +  F3  a3  .  . 

v3i  (0  ,  0.25)  £33(Oc  0.25,  t)  =  (Fi  -  Gi)  a& 


+  (F2  -  G3)  a2  +  (F3  -  GO  a3 

a  ci 


v31  (0  ,  0.50)  e33  (0  ,  0.50,  t)  =  (Fi  -  Gi)a. 


(7) 
(8) 
(9) 

(10) 


+  (F2  -  G3)  a2  +  (F3  -  GO  a, 


(11) 


v31  (0  ,  0.75)  £33  (0  ,  0.75,  t)  =  (Fi  -  Gi)  a 
c  c  t. 


+  (F2  -  G3)  a2  +  (F3  -  Git)  a3 


(12) 


Fi,  F2  and  F3  can  be  found  by  Cramer's  rule  or  other  methods  from  Eqs. 
7   to  9,  whereas,  Gi,  G3  and  G^  can  be  solved  from  Eqs.  10  to  12.   Sub- 
scripts a,  b,  and  c  correspond  to  stress  levels  of  0.25,  0.50  and  0.75 


respectively.   Once  these  6  kernels  are  found,  the  complete  time  history 
of  the  strains  are  (8) 

e33  (0c,  S,  t)  =  J*  Fx  a    (d)  dd  +  /*/*  F2  a  (Ci)  a  (?2)  d^  dc2 

+  Illlll   F3  a  (d)  a  (C2)  a  (?3)  d?j  d^2  dC3     (13) 

e22  (0c,  S,  t)  =  en(0c,  S,  t)  =  /*  (Fj  -  Gj)  a  (d)  d^ 

+  /0/o  [(F2  "  G3)  °  (Cl)  °  (C2)]  dCl  d^2 

+  /oOo  [(F3  "  Gt+)  '°    (Cl)  °  (?2)  °  (^)] 

x  dCi  d^2  d^3 (14) 

e32  =  £31  =  £21  =  0    (15) 

£l>  ^2>  ^3  are  dummy  variables. 


EXPERIMENTAL  PROCEDURE 


Sample  Preparation 


One  of  the  main  difficulties  in  mechanical  testing  of  an  oil  shale 
is  the  preparation  of  samples.   Conventionally,  all  the  oil  shale  sam- 
ples were  cut  mechanically  by  diamond  saw  blades  and  cores.   To  avoid 
any  micro-cracks  in  the  samples  produced  by  cutting,  a  wire  saw  with 
diamond  impregnated  wire  was  used.   The  wire  saw  cuts  smoothly  and 
accurately,  with  a  minimum  loss  of  kerf. 

Since  oil  shale  is  a  stratified  rock,  care  was  taken  to  select  as 
uniform  samples  as  possible.   Layers  of  obviously  foreign  composition, 
major  elastic  inclusions,  and  locations  where  faults  existed  were  avoided, 
The  set  of  3  to  4  duplicate  oil  shale  samples  required  for  creep  test- 
ing, obtained  by  exploiting  the  lateral  uniformity  of  oil  shale,  were 
cut  from  a  single  horizontal  (3.175  cm.  thick)  layer  from  a  core  of  the 
upper  Green  River  Formation.   The  test  sets  were  selected  with  differ- 
ent oil  yields.   The  size  of  each  specimen  used  for  creep  testing  was 
1.27  cm.   square  in  cross-sectional  area  and  3.175  cm.  long. 

Instrumentation 

Prior  to  creep  testing,  a  Tinius  Olsen  universal  testing  machine 
was  used  to  determine  the  ultimate  compressive  strength  of  the  material 
of  the  same  horizon  as  the  materials  to  be  tested  for  creep. 

For  creep  tests,  specimens  were  placed  between  spherical  seat  and 
the  fixed  plate  of  the  creep  frame.   The  creep  machine,  (Fig.  3)  consists 
of  a  Conbel  Model  355  testing  machine,  pneumatic  driven.   Dial  gages 
readable  to  0.00254  mm  (.0001")  were  used  to  measure  the  displacement  of 
the  specimens. 


Figure  3.   Creep  Machine 


Three  to  four  duplicate  samples  were  tested  with  various  loads, 
ranging  from  25%  through  75%  of  the  ultimate  compressive  strength  of 
the  material  in  the  layer.   In  this  manner,  the  stress  effect  can  be 
determined  on  the  specimens  of  the  same  oil  and  mineral  content. 
Specimens  were  placed  between  spherical  seat  and  fixed  plate  of  the 
creep  frame.   The  spherical  seat  resting  on  a  lubricated  socket  was 
used  to  insure  a  uniform  contact  pressure  during  testing.   An  initial 
load  of  approximately  890  N  (200  pounds)  was  applied  to  check  for  eccen- 
tricity and  for  specimen  alignment.   The  exact  load  applied  to  the 
specimen  was  read  from  the  dial  gage  of  the  creep  machine  which  was 
calibrated  before  it  was  used  as  a  loading  device.   Also,  a  dial  gage 
readable  to  0.00254  mm  (0.0001")  was  placed  between  a  fixed  bar  and  the 
spherical  seat,  and  the  air  was  introduced  slightly,  just  enough  to  set 
the  dial  gage  into  motion.   At  this  instant,  the  exact  dial  reading  was 
taken.   The  load  was  then  increased  to  full  load  and  maintained  through- 
out the  testing  by  an  air  pump. 


Prior  to  creep  testing,  the  specimens  were  surface-dried,  weighed, 
measured,  and  then  placed  in  the  creep  frame  for  testing.   The  vertical 
displacement  of  the  specimen  was  measured  after  the  initial  load  was 
applied  and  also  when  the  full  load  was  applied.   The  time  interval 


8^ 

for  measuring  the  displacement  depended  on  the  displacement  rate  of  the 
specimen  tested. 

Due  to  the  small  displacements,  great  precaution  was  taken  in  meas- 
uring them.   Typical  data  of  creep  strain  versus  time  is  shown  in  Fig.  1, 
It  shows  the  creep  behavior  with  an  initial  non-linear  range  and  gradu- 
ally assumes  a  constant  strain  rate. 


TEST  RESULTS  AND  THEIR  STATISTICAL  ANALYSIS 

Data 

The  following  variables  were  found  significant  in  the  statistical 
analysis: 

0   =  organic  content  in  fraction  by  volume;  and 

S   =  stress  level  as  a  percentage  of  ultimate  compressive  stress. 

From  the  oil  yield  values  M  in  gallons  per  ton  (1  gallon/ton  = 
4.19  cm3/kg)  determined  by  Fischer  Assay,  the  organic  content  in  fraction 
by  volume,  0  ,  was  calculated  using  the  relation  derived  by  Smtih  (10): 

n         157. 76M  nfi. 

c  "  (0.957M  +  107)  100 ^      ' 

Statistical  Analysis 

From  the  raw  data,  it  was  found  that   n  varied  linearly  with  re- 
spect to  0C  and  S  (Figs.  4,  5);  that  e+  varied  linearly  with  0C  and 
parabolically  with   S  (Figs.  6,  7);  that  e°  varied  parabolically  with 
0  but  linearly  with  S  (Figs.  8,  9).   Thus  the  form  of  the  regression- 
al  equations  was  determined. 

Based  on  regressional  analysis,  the  following  equations  were  ob- 
tained which  relate  the  mechanical  properties   n,  e+,  and  e°  ,from  Eq.  (1) 
with  the  organic  contents  and  stress  levels: 

n  =  0.2801  -  0.8245  0  -  0.1750  S  +  1.5827  0  S (17) 

c  c 

Coefficient  of  determination,  r2  =  0.7440 
Number  of  samples,  N   =  13 

Standard  error  of  estimate,    s   =  0.0286 


Mean  values,  n   =  0.1736 

0C  =  0.1870 
S  =  0.4808 


Fig.  4.   Parameter,  n  vs  Organic  Volume,  0, 
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e+  =  -0.0005  +  0.0027S  +  0.0004  Oc  -  0.0026S2  +  0.0034  OcS   (18) 

r2  =  0.8383 

N  =  13 

s  -  0.0001 

7+  =  0.0004 
Oc  =  0.1612 
S  =  0.4808 

e°  =  0.0141  +  0.0885  Oc  -  0.4130  02  +  0.1085  0CS (19) 

r2  =  0.8462 

N  =  12 

s  =  0.0023 


e_°  =  0.0225 
Oc  =  0.1627 
S  =  0.4808 


Typical  prediction  lines  based  on  Eqs.  17  to  19  are  plotted  in 
Figs.  4  to  9  as  solid  lines. 


CONCLUSIONS 

The  mechanical  properties  of  oil  shale  under  creep  tests  were  in- 
fluenced by  the  volumetric  organic  content  and  stress  levels.   The  or- 
ganic content  is  an  important  parameter.   The  cross  product  of  organic 
content  and  stress  level  also  improved  the  multiple  regression  fit.   The 
creep  can  be  represented  by  a  nonlinear-strain-time  model  as  indicated 
by  Eq.  1,  yielding  the  parameters  e°,  e+  and  n.   From  these  parameters, 
a  nonlinear  pseudo  3-D  creep  model  can  be  built,  relating  the  strains  to 
the  nonlinear  stresses. 
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